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ABSTRACT 

We use a grid-based shallow water model to simulate the atmospheric dynamics of the transiting hot 
Jupiter HD 209458b. Under the usual assumption that the planet is in synchronous rotation with zero 
obliquity, a steady state is reached with a well-localized cold spot centered 76° east of the antistellar 
point. This represents a departure from predictions made by previous simulations in the literature 
that used the shallow water formalism; we find that the disagreement is explained by the factor of 30 
shorter radiative timcscalc used in our model. We also examine the case that the planet is in Cassini 
state 2, in which the expected obliquity is ~ 90°. Under these circumstances, a periodic equilibrium is 
reached, with the temperature slightly leading the solar forcing. Using these temperature distributions, 
we calculate disk- integrated bolometric infrared light curves from the planet. The light curves for the 
two models are surprisingly similar, despite large differences in temperature patterns in the two cases. 
In the zero-obliquity case, the intensity at the minimum is 66% of the maximum intensity, with the 
minimum occuring 72° ahead of transit. In the high-obliquity case, the minimum occurs 54° ahead of 
transit, with an intensity of 58% of the maximum. 

Subject headings: stars: planetary systems - hydrodynamics 



1. INTRODUCTION 

Hot Jupiters orbit their par ent stars with periods P < 7 
days (Mayor & Queloz 19951. They therefore receive a 
large stellar flux and present an opportunity to study plan- 
etary behavior under conditions not seen in the solar sys- 
tem. Furthermore, these planets are expected to have 
evolved tidally to spin-synchronicity. In the event that 
the planetary spin and orbital angular momentum vectors 
are aligned, the planet keeps one hemisphere in perpet- 
ual daylight and the other in perpetual darkness. The 
combination of intense irradiation and a permanent night 
side provides interesting modelling challenges t hat can be 
tested with observations in the infrared band (Charbon- 
neau et aL] [2005; Do ming et al.|[2005] |Harrington et al.' 
2006P - 

The transiting extrasol ar 
(Charbonneau et al.||2000 



giant planet 
Henry et al.|[2000h 

tunity to study th 



HD 209458b 
provides a 

particularly suitable opportunity to study these unique 
dynamics. It orbits its primary at a distance of 0.045 
AU with a period of 3.52474541 days (Wittenmycr et al 
2005}. The planet h as a mass of 0.697^7 
oTT32i?j ± 0.03i?j (iCharbonneau et al. " 



anc 



a radius 
20061, ~ 30% 



lar ger th an current planet ary models predict (Eaughlin et 
al. 2005 1 . As reviewed by Charbonneau et al. ( 2006 ), the 
planet must have a large source of interior heat; several 
explanations have been advanced to explain how this ex- 
tra heat is produced. In this Letter, we investigate the 
potential observat i onal c onsequences of the hypothesis of 
Winn & Holman| (|2005|) 



who suggest that HD 209458b 
assini state 2, in which the orbital 



may be trapped in 
angular momentum vector and the planetary spin vector 
precess at the same rate and the obliquity is forced to ap- 
proach 90° during planetary migration. In this scenario, 
the anomalous radius results from the heating caused by 
ongoing obliquity tides. In this case, the atmospheric ir- 
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radiation pattern is very different from that of the zero- 
obliquity case. It is therefore interesting to investigate 
whether this difference results in a qualitatively different 
infrared light curve for the planet. 

Several hydrodynamic models have been used to inves- 
tigate the atmosphere of HD 209458b under the assump- 



tion that the planet has zero obliquity (Cho et al 



Cooper fc Showman 2005 
fc ShowmanT f2005) r 
dimensional global ci 



mo 



de 



C 



2003 



ooper 



Burkert et al.||2005[) . ^ 

the atmosphere with a three- 
global circulation model. They find that at 
atmospheric optical depth r ~ 1, a superrotating, super- 
sonic equatorial jet advects the hottest region of the at 



Cho et al. (20031 



mosphere ~ 60° downwind. In contrast, 
use a shallow water model with a radiative time constant 
Trad — 10 days. They obtain circumpolar vortices that 
sequester cold air and maintain a local temperature vari- 



ati on of up to 1000 K. However, subsequent work (Iro et 
aL||2005| ) has suggested that the photospheric r ra d should 
be on the order of hours, rather than days. At present, 
there seems to be little agreement between models. For- 
tunately, because the atmosphere is expected to exhibit 
extreme temperature variations, it is possible that obser- 
vations of the planet's infrared signature will be able to re- 
solve these discrepancies. Therefore, a second goal in this 
Letter is to investigate t he cause of t he diff eren ce between 



Showman (20051. 



our results and those of Cho et al. (20031 and Cooper & 



2. ATMOSPHERIC MODEL 

Our simulations of the atmospheric flows are computed 
with a sh allow water model similiar to that used in lCho etl 
al. (2003). This model portrays the atmosphere as a thin 
layer of incompressible fluid at a constant p and therefore 
is a reasonable approximation for subsonic stratified flows. 
Framed in terms of the temperature T and the horizontal 
wind velocity v = ucf> + v9, the shallow water equations on 
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an irradiated, rotating sphere are 
dT 
~dt 
<9v 

at 

where R is the specific gas constant, / = 2wsin6* is the 
Coriolis parameter, h is the surface normal unit vector, 
and Ft is the thermal forcing. A rigorous treatment of the 
thermal forcing would involve a radiative transfer scheme; 
in this initial investigation, we choose instead to follow 
Cho's parametrization of the heating. We assume an equi- 
librium temperature distribution 



AT 

T eq = T ef f + — cos a, 



(3) 



where T e ff is the equilibrium blackbody temperature, 
AT is the maximum equilibrium day/night temperature 
difference at the layer being modeled in the absence of 
any hydrodynamic effects, and a is the angle between a 
given point (0, 9) on the planet's surface and the substel- 
lar point. The equilibrium day/night temperature differ- 
ence is a free parameter, and corresponds to the efficiency 
with which the planet can redistribute heat at an adia- 
batic depth. This AT parameter should not be confused 
with the actual temperature variations that develop as 
the atmosphere evolves. Rather, it measures the magni- 
tude of the difference in forcing between the day and night 
hemispheres. For the zero-obliquity case, we simply have 
cos a = cos 4>cos 9. We then impose Newtonian relaxation 
to this equilibrium temperature, so that 



Fx 



T — T e , 

7~rad 



(4) 

= 8 hour s, the ex- 
200 mbar (Iro et al. 



We adopt a radiative timescale r ra 
pecte d value at a pressure depth p = 
2005 1. For the results presented here, we take T e // = 1400 
K, appropriate to a Bond albedo A ~ 0.2. We further as- 
sume AT = 800 K. This relatively large value is motivated 
by the recent full-phase ob servations of v Andromedae b 
by |Harrington et al ( 2006 ) . 

The shallow water equations are integrated using a vec- 
tor spherical harmonic t r ansfor m method proposed by 
|Adams fc Swarztrauber ( |1999[ ). The models are run 
at resolutions between 256 x 127 (T85) and 512 x 257 
(T171). Numerical stability is mai ntained using fourth- 



order hyperviscosity ( Fornberg|1998 1 , which pre vents spec- 
tral blocking at high wavenumbers ( |Boyd|[200l| . We start 
with an isothermal (T = T e f f) atmosphere initially at rest, 
and integrate until either a steady state or a periodic equi- 
librium is reached, a process that is complete within 250 
rotation periods at the radiative timescales considered. 

3. RESULTS OF THE ZERO-OBLIQUITY CASE 

In the case that the planet has negligible obliquity, our 
model predicts that the atmosphere will exhibit an ap- 
proximately steady state condition, as seen in Fig. T] The 
dominant feature is a localized cold spot on the night side. 
However, this cold spot is not directly at the antistellar 
point; rather, its center is 76° to the east; the coldest point 
is 63° east of the antistellar point. Interestingly, there ap- 
pears to be a strong westward equatorial jet, with an av- 
erage zonal wind speed of 780 m s _1 at the equator. The 



eastward displacement of the cold spot then results from 
the shallow water analog of the Bernoulli effect, rather 
than simple downwind advection of cold air. The temper- 
ature distribution is quite asymmetric, with the maximum 
110 K warmer than T e ff and the minimum 280 K colder 
than T e ff 



This pictu re differs dramatically from the results of Cho 
|et al.| (|2003 h , despite the similiarity of the model. Cho et 
|al.| ( |2003p start with random initial vorticity and use a long 
radiative time constant r ra d = 10 days. After some tens 
of rotation periods, a dynamic equilibrium has evolved in 
which circumpolar vortices produce cold spots that orbit 
about the poles. This leads to a situation where the tem- 
perature minimum can, at times, be located on the day 
side of the planet. 

While we start the atmosphere at rest, rather than using 
vortical initial conditions, the primary differ ence between 
our simulation and that of Cho et al. (20031 is the radia- 
tive time constant; a value of T ra d = 8 hr is now believed 



to be more realistic (Iro et al. 20051. Indeed, their r ra d 
is sufficiently long that the forcing appears to have only 
negligible impact on the evolution of the atmosphere. We 
ran simulatio ns starting f rom v ortical initial conditions as 
described by Cho et al. (20031 with no thermal forcing; 
within 40 rotation periods, the flow pattern was virtually 
identical to the Cho result with forcing. The end state of 
these simulations may be seen i n Fig. |2| The con clusion is 
that the atmospheric flow that Cho et al. ( 2003 1 describe 
is more influenced by their initial conditions than it is by 
the intense stellar radiation. 

While our result s are more similar to those of |Cooper| 



ancies 
mension 



& Showman (2005 1, there ar e still some notable discrep- 
Cooper & Showman (20051 use a fully three di- 



global circulation model, starting with a cold 
isothermal atmosphere at rest; after some time, a steady 
state emerges in the upper atmosphere. They find that 
a strong eastward jet develops at the equator, advect- 
ing the hottest region ~ 60° downwind. In contrast, 
we have a strong westward jet at the equator, with the 
cold spot being pulled upwind. The reason for this differ- 
ence is not immediately clear. There is some indication 
that the development of superrotating jets require s three- 
dimensional processes ( Cooper fc Showman||2005[ ) . If this 
is the case, they clearly would be unable to arise in our two- 
dimensional shallow water model. On the other hand, a 
two-dimensional model allows higher horizontal resolution, 
which could reveal processes that would not be resolved in 
a full three-dimensional model. Since the atmospher e at 
the infrare d photosphere is expected to be radiative (Iro 
et al.|2005 1 , vertical mixing should be of secondary impor- 



tance compared to the horizontal flow, and the use of a 
two-dimensional model is justified. 

4. RESULTS OF THE HIGH-OBLIQUITY CASE 



Winn & Holman (120051 suggest that the large radius 
of HD 209458b could be explained if the planet were in 
Cassini state 2. A Cassini state is a spin-orbit resonance 
in which the planetary spin axis and the orbit normal pre- 
cess about the same axis at the same rate. In Cassini state 
2, they are on opposite sides of the precession axis, and the 
planet must have non-zero obliquity. It is expected that 
as the planet migrates inward and spin-synchronization 
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occu rs, the axial tilt will in crease until it is very close to 
90° ( [Winn fc Hoiman|[2005l ) . 

In this case, the entire surface of the planet will receive 
stellar radiation; the time-dependent stellar elevation fol- 
lows 

cos a — cos 9 cos 9* (cos </> cos </>* + sin <f> sin </>* ) +sin 9 sin 9* , 

(5) 

where = — cut and 9* = —tot give the position of the sub- 
stellar point on the planetary surface. The forcing itself is 
still given by equations [3] and [4] As may be expected, the 
atmospheric dynamics in this case differ greatly from those 
in the zero-obliquity case. Here, a periodic equilibrium is 
reached, with one period per rotation. The atmospheric re- 
sponse leads the forcing by ~ 55° , so that the temperature 
begins to rise 13 hr before the onset of direct solar heating. 
Temperature maxima are higher at the poles than at the 
equator as equatorial winds redistribute the heat. Also 
notable are colder spots on the illuminated face produced 
by ring-shaped vortices that form at the poles and then 
break as they move to lower latitudes. The temperature 
distribution is considerably more symmetric than that of 
the zero-obliquity case; the maximum and the minimum 
are ~ 250K warmer and colder than T e //, respectively. 

5. INFRARED SIGNATURES 

Having determined the photospheric temperature, it is 
possible to calculate the relative intensity of the planet's 
infrared emissions. We assume that the planet radiates as 
a blackbody; then the intensity of radiation emitted is 



(n-f)fcT 4 cos 9A9A0, 



(6) 



where h is a unit surface normal vector, r is a unit vector 
directed from the planet towards earth, and A: is a nor- 
malization constant chosen so that I = 1 at maximum 
intensity. We sum over all grid points on the portion of 
the planet facing earth (i.e., where h ■ r > 0). 

The resulting light curves for both zero obliquity and 
high obliquity are seen in Fig. [4] Due to the temperature 
asymmetry apparent in the zero-obliquity case, the max- 
ima are considerably wider than the minima. This is the 
major contrast between the two cases; in the case of high 
obliquity, the maxima and minima are symmetric. How- 
ever, neither curve diverges significantly from a sinusoidal 
shape. In both cases, the extrema lead the forcing. In 
the zero-obliquity case, the minimum has a phase offset of 
72°, while the maximum is offset by only 43°. The high- 
obliquity case is closer to symmetric, with slightly smaller 
phase offsets of 54° and 58° at the minimum and maxi- 
mum, respectively. 

These light curves are difficult to reconcile with the re- 
cently repo rted data from v Andromedae b (Harrington 
et al.|2006 l, in which the phase offset between transit and 
the infrared minimum is much smaller and can be inter- 
preted as arising from a situation in which very little heat 
is transferred from t he day side to the night side. Har- 



rington et al. ( 2006 1 obtain a best fit with a phase shift 
of 11°, and the data are not inconsistent with zero phase 
shift. We suggest two possible causes of this small phase 
shift. First, at the 24 /im photosphere from which the data 
were taken, the pressure depth is quite low due to infrared 
absorption by water molecules. Indeed, the 24 -/xm photo- 
sphere occurs at a pressure depth of 20 mbar (Fortney et 



aL] |2005[ ), at which r rad = 5.7 hr fllro et al.|[2005| . Our 
model predicts a smaller phase offset as T ra d decreases. 

In addition, our model indicates that lower atmospheric 
temperatures will also lead to smaller phase offsets. Using 
planetary parameters appropriate to v Andromedae b, an 
atmospheric temperature T e ff — 1500 K and a radiative 
timescalc of 4 hr yields a phase offset of 49°; decreasing 
the atmospheric temperature to 800 K reduces the phase 
offset to 23° . While more careful determination of the pa- 
rameters and a more rigorous treatment of the radiative 
forcing will be necessary to obtain a good fit to the v An- 
dromedae b data, it seems quite plausible that the small 
phase offset is the result of a combination of the smaller ra- 
diative timescales and lower temperatures associated with 
the low pressure depths at which the measurements were 
made. 

6. CONCLUSIONS 

We have presented the results of a two-dimensional grid- 
based hydrodynamic simulation of the atmosphere of the 
transiting extrasolar giant planet HD 209458b. In addi- 
tion to the traditional case in which the planet's axial tilt 
is close to zero and the same side always faces the star, we 
have also considered the case that HD 209458b is tilted 
on its side, with an obliquity near 90°. This configura- 
tion is dynamically stable, and is motivated by the need 
to explain the planet's anomalously large radius. 

In both cases, a stable equilibrium develops, although 
the specific temperature distribution is markedly differ- 
ent. In the zero-obliquity steady state temperature 
profile is reached, characterized by a diffuse hot region 
blown around the globe by eastward jets at higher lat- 
itudes and a localized cold spot produced by a Bernoulli 
effect within an overall westward jet at the equator. In the 
high-obliquity case, the temperature distribution is dom- 
inated by the contribution of radiative forcing, although 
the temperature extrema lead the forcing by ~ 13 hr. 

In both cases, we expect to see observable changes in the 
planet's infrared spectrum; the maximum infrared emis- 
sion is between 50% (zero obliquity) and 70% (high obliq- 
uity) greater than the minimum . However, the difference 
in the light curves between the two states is likely to be 
difficult to observe experimentally, as both light curves are 
approximately sinusoidal. In addition, the phase lag be- 
tween illumination and infrared emission is similiar in both 
cases. 

While the results of this preliminary investigation do not 
provide a clear observational distinction between the low- 
and h igh-obliquity cases, the results of |Harrington et al~ 
(2006| demonstrate that analysis of the infrared emissions 
of extrasolar planets is a present possibility. We believe 
that improvements to our model currently being tested 
will enable us to provide a good fit to the Harrington data 
and to other forthcoming observations of the infrared emis- 
sions of extrasolar planets. 

We acknowledge support from the NASA Plane- 
tary Geology and Geophysics program through grant 
NNG04GK19G, and from the NSF through CAREER 
grant AST-0449986. The code used for the simulations 
included the SPHEREPACK 3.0 subroutines provided by 
the University Corporation for Atmospheric Research. 
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Fig. 1. — Steady state temperature of HD 209458b in the zero-obliquity case. Day side temperature distribution is shown on the left; night 
side on the right. 
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Fig. 2. — Dynamic equi librium reached fr om the evolution of vortical initial conditions in the absence of forcing. Of note is the similarity 
to the results presented in |Cho et al.| |2C)03| |. 
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Fig. 3. — Time-dependent temperature of HD 209458b in the high-obliquity case. Each successive frame advances by one-fourth of a 
rotation period; the temperature variation is periodic, repeating once per rotation. The top row shows an equatorial view, while a polar view 
can be seen on the bottom row. 
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Fig. 4. — Infrared emissions from HD 209458b over one full rotation period. Transit occurs at t = 0. The solid line shows the light curve 
from the zero-ob liquity case; the dashed line shows the light curve from the high-obliquity case. The light curve obtained by |Cooper fc| 
|Showman| ( 2005 ) , shown by the dotted line, is included for comparison. 



